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Carnegie  Mellon  University.  Pittsburgh.  PA 


i Received  October  II.  1991) 


Abstract  The  anisotropic  third-order  nonlinear  optical 
susceptibility  of  monodotnain  samples  of  nematic  sulutions 

of  poly{  1 ,4-phenylene-2,6-benzobisthiazole  1,  PBT.  a  rodlike 
polymer,  have  been  studied  by  third  harmonic  generation.  The 
principal  component  of  x(3)(-3co;a>,o>.w)  is  not  along  the  nematic 
director,  as  opposed  to  the  linear  dielectric,  which  has  the 
maximum  refractive  index  along  the  director. 


INTRODUCTION 

Nonlinear  optical  properties  of  organic  material?  have  been  of 
intense  interest  for  the  past  decade  and  more,  e.g.,  see  references 
1-8.  The  susceptibilities  to i are  macroscopic 

properties;  their  relation  to  molecular  characteristics  is  discussed 
below.  In  the  dipole  approximation,  the  macroscopic  polarization 
P(to)  of  a  material  under  the  influence  of  an  applied  fields  Efwi), 
E(co2),  etc.,  propagating  with  frequencies  to],  0)2,  etc.  (e.g..  optical 
radiation)  may  be  expressed  in  Cartesian  coordinates  as;1'** 

P.(oo)  =  eofx/j”  Ej  +  xsjk  EjEk  +  XijS  EjEkEi  +  ...]  ( 1 ) 

where  co  =  coj  +  0)2+  003,  etc.;  the  component  notation  in  Eqn.  (I),  in 
which  the  dependence  on  frequency  is  suppressed  for  convenience, 
will  be  used  below.  In  this  study,  we  report  measurements  of 

9  3  *4  0  1  2  8  8 
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3o>:co. co. 0) )  by  third  harmonic  generation  on  a  fully  aligned 
nematic  solution  of  the  rodlike  polymer  poly[  1 ,4-phenvlene-2.6- 
benzobisthiazoiel.  PBT 


PBT 


The  poivmer  is  dissolved  in  a  strong  protic  acid  and  is. 
therefore,  protonated.u  Prior  measurements  of  y1  ^Voutco.oj.oj )  have 
been  for  isotropic  solutions  of  PBT  in  the  same  solvent  and  on  PBT 
dispersed  in  an  organic  matrix.  ‘  1 1  In  addition.  xl  3  ’< -Jtoico.to.u)  i  and 

-or.ta.-oj.to )  have  been  studied  for  undiluted  PBT  i  x'  '( -ar.co.-ui.ci ) 
being  derived  from  degenerate  lour  wave  mixing;.11-1-  The  aligned 
solutions  studied  iiete  afford  the  opportunity  to  elucidate  certain  of 
the  components  lx’  '^ijkl  of  the  susceptibility  tensor.  The  PBT 
samples  studied  previously  did  not  possess  the  alignment  of  the 
samples  examined  in  this  works.  For  an  isotropic  sample,  third- 
harmonic  measurements  give  the  averaged  value  fx(3)hso  = 
rix<3>iuij/5  o 

Previous  studies  with  the  present  nematic  PBT  solutions  have 
shown  them  to  be  strongly  birefringent,  with  a  birefringence  An 
that  increases  essentially  linearly  with  the  the  polymer  volume 
fraction.14  This  property  was  successfully  exploited  to  extract  a 
value  for  the  order  parameter:  S  =  0.94.  Such  a  value  reflects  the 
very  high  alignment  of  the  polymer  rods  in  the  present  solutions.14 
The  refractive  index  measurements  show  that  both  the 
extraordinary  and  the  ordinary  refractive  indices,  ng  and  no, 
respectively,  depend  on  the  wavelength  X,  with  (X  is  E  or  O) 


( nx~  +  2 ) 


(2) 


where  ae  and  Ao  were  found  to  be  250  and  200  nm.  respectively, 
and  K  was  0.27  for  both  ng  and  no-14  Analysis  of  the  experimental 
results  obtained  here  makes  use  of  these  refractive  indices. 
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n:9)'4< 


EXPERIMENTAL 

Materials:  The  PBT  (Mw  =  34,000)  was  provided  by  SRI 

International.  The  poivmer  was  dried  in  vacuum  and  dissolved  in 
distilled,  anhydrous  methane  sulfonic  acid  (MSA)  using  procedures 
discussed  elsewhere.9  Dissolution  in  MSA  is  accompanied  by 
protonation  of  the  solute  macromolecules.9  Solutions  of  PBT  are 
nematic  when  formed  with  a  volume  fraction  o  above  a  critical 
value  <pni.  where  o  \ \  depends  on  the  chain  length  and  the 

temperature:  <pni  =  0.03  for  the  PBT  used  at  25°C.  The  alignment 
procedure,  described  elsewhere.14-15  consists  of  a  surface  alignment 
preparation  using  a  suitable  How  in  a  rectangular  channel,  followed 
by  imposition  of  an  external  magnetic  field  (5  to  ''  Tesla  >  in  the 
sample  plane,  and  along  the  flow  direction,  to  ^peeu  the  bulk 
alignment.  With  suitable  surface  alignment,  the  nematic 
monodomains  are  stabilized  for  an  indefinite  duration. 

Nonlinear  Optical  Measurements:  Third  harmonic  generation 

(THG)  technique  has  been  used  to  determine  the  third-order 
susceptibility.1-2-11  The  THG  apparatus  utilizes  a  Raman  Cell  filled 
with  methane  gas  to  provide  a  fundamental  intensity  I,0at 
wavelength  X  ~  1542  nm,  when  pumped  at  X  =  1064  nm  by  a 
pulsed  Nd:YAG  laser,  generating  third  harmonic  intensity  at  X  = 
514  nm.  The  Maker  Fringe  Pattern  (MFP),  generated  by  h(1)  and 

used  to  study  the  third-order  susceptibility,  may  be  formed  on 

rotation  of  a  plane  parallel  sample,  or  translation  of  a  wedge- 
shaped  sample,  in  an  incident  beam  1^, 2,3.7  The  incident  beam  is 
split  into  two  beams,  directed  to  the  sample  and  the  reference  to 
account  for  fluctuations  in  the  incident  intensity.  A  weakly  focused 
beam  generates  a  sufficiently  intense  signal  for  accurate  detection 

with  the  samples  of  interest  here.  A  Fresnel  Rhomb  in  the  incident 
beam  permits  control  of  the  incident  polarization,  and  an  analyzer 
in  the  THG  beam  allows  selection  of  the  polarization  of  the  output 
signal. 

bxpeiimcms  used  plane  parallel  ceils,  provided  by  Helima  Cells 
Inc.  The  cell  walls  were  specified  to  X/4  surface  flatness:  the 
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sample  thickness  L  was  375  urn.  The  reduced  intensity  R(8.L)  is 
given  by 


R(e.u 


Ij^Q.DsAMP 

r3U(0X)REFR 


(3) 


where  SAMP  and  REFR  <tand  t'or  the  sample  and  reference, 
respectively,  and  8  is  the  incidence  angle  and  E  the  physical 
thickness;  8  is  varied  by  rotation  of  the  sample  t  fixed  Li  in  the  use 
of  a  parallel  slab,  and  L  is  \aned  by  translation  of  the  sample  (fixed 
8,  essentially  equal  to  zero)  in  the  use  of  a  wedge-shaped  sample.2 
A  fused  silica  plate  <L  =  1.00  mm)  was  used  as  the  reference. 

The  MFP  exhibits  damped  and  symmetric  oscillations  about  8  = 
0  for  the  rotating  parallel  slab  cell.2b-c  A  general  expression  for 
R(0.L)  for  an  isotropic  material  is  given  by:'"1 


R(0,L)  =  Kr 


1Q(8)W(0.L)1samp  r 
[Q(0)W(0,L)1refr  j 


(4) 


W(0.L) 


IX 


(3), 


In¬ 


to 


FFF 


Qie.L) 


(5) 


where  Ijc<3)Ieff  is  the  effective  component  of  xl3).  e*g-,  see  Eqn  (1), 
Q(9)  is  a  weakly  decreasing  function  of  increasing  angle,  related  in 
part  to  the  Fresnel  factors  (Q(0)  is  unity),2  Kr  is  the  ratio  between 
the  incident  beam  directed  to  the  sample  and  that  directed  to  the 

reference,  and  n^  and  n^  are  the  refractive  indices  at  the 

fundamental  and  the  THG  frequencies,  respectively.1  4 

The  function  Q(0,L)  accounts  for  the  interference  effects 
between  two  waves,  propagating  with  frequency  3<u  and  velocities 
c/n3oj  and  c/n^.  giving  rise  to  the  MFP.  For  the  case  with  weak 

absorption,  characterized  bv  absorbances  cto,  and  for  the 

incident  and  THG  beams,  respectively,  Q(0,L)  may  be  expressed  as:2 


fl2(e,L)  =  2b 


2b 


sin:l\g(0,L)] 


( 6a) 
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Here, 

V(6.L)  =  jiL/2Lc(0)  (6b  I 

4a  =  exp(-a3co)  +  exp(-2a(„)  (6c) 

2b  =  exp(-a3{0/2)exp(-a(l))  (6d) 


where  a  =  |iL/cos0’,  with  p  the  extinction  coefficient,  and  the 
coherence  length  Lc(0)  may  be  expressed  as 

Lc(0)  =  X/6ln3Wcos<  0\U))-n„)cos(0'(J))i  (1) 

with  Lc(0)  =  X. / 6 1  n 3 o> - n a, ! ;  the  primes  denote  propagation  angles 
measured  in  the  sample,  related  to  angles  in  air  through  SnellT  law. 
using  the  appropriate  refractive  index,  see  below.  The  effect  of 
weak  focusing  of  the  incident  beam  is  neglected  as  too  smail  to  be 
important  under  the  conditions  used.2  In  practice,  the  parameter 
Kr  is  evaluated  by  determination  of  Rstd(O)  for  the  response  to  a 
standard  with  known  I*  Icff.  n3w.  Rw,  and  thickness  L.  A  BK7  plate 
(L  =  3.00  mm)  was  used  as  the  standard  in  the  cited  study.  It 
should  be  noted  that,  for  a  finite  Lc(0),  R(0,L)  may  take  on  any 
value  from  zero  to  its  value  Rmax.  as  \|/(9,L)  is  not  usually  equal  to 
an  odd  multiple  of  Jt/2  at  0=0.  However,  siny  may  be  taken  as  unity 
for  the  successive  maxima  in  the  fringe  pattern,  providing  a 
method  to  estimate  of  lx(3\FF/'n3o>“  -  nu2l  from  data  on  R(0,L).  We 
have  also  shown,  in  previous  measurements,  that  the  THG  analyzed 
here  is  mainly  generated  from  the  solute  polymer.  The 
contribution  from  the  solvent  and  the  cell  walls  are  negligible:  less 
than  5%  for  solutions  with  volume  fraction  ip  =  0.02  -  0.055. 10 

For  an  anisotropic  material,  the  parameter  lx  1i;ff  becomes  the 
modulus  of  an  appropriate  component  of  the  tensor  x(V,  depending 
on  the  polarization  of  the  light  and  the  orientation  of  the  director  n. 
In  addition,  the  birefringence  expected  with  an  anisotropic  material 
may  modify  the  MFP,  see  below. 

i 
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FIGURE  I  Maker  Fringe  Patterns  tor  a  nematic  solution  of  PBT. 
with  parallel  polars  and  the  incident  polarization  orthogonal  to 
the  director  (e.g.,  arrangement  V-h-V  or  H-v-H). 


RESULTS 

With  nematic  solutions,  the  response  R(8,L)  depends  on  the 
polarizations  of  the  incident  light  and  the  detected  THG  signal,  and 
the  orientation  of  these  relative  to  the  sample  director  (the 
preferred  optical  axis).  Several  geometric  arrangements  are  of 
interest,  depending  on  the  relative  orientations  of  the  polarizations 
of  the  incident  and  THG  fields,  Einc  and  Ethg.  respectively,  and  the 

director  n  of  the  monodomain  with  respect  to  the  axis  of  rotation; 

the  latter  was  always  vertical  in  the  arrangements  used.  For 

convenience,  the  dependence  of  R(0,L)  on  L  is  suppressed  hereafter. 
All  data  were  obtained  using  1542  nm  incident  wavelength.  The 

large  birefringence  of  the  PBT  monodomain  complicates  the 
analysis  is  some  cases,  as  the  relevant  refractive  index  varies  with 
0  for  certain  arrangements,  see  below.  In  the  following,  the 
notation  gives  the  orientations  of  these  fields  in  the  order 
ETHG;n;EfNc<  e.g.,  V-h-H  signifies  an  arrangement  with  Ethg  vertical 
(parallel  to  the  rotation  axis)  and  both  Einc  and  n  horizontal. 
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FIGURE  2  Maker  Fringe  Patterns  for  a  nematic  solution  of  PBT,  with 
parallel  polars  and  the  incident  polarization  parallel  to  the 
director  arrangement  V-v-V  (a)  and  H-h-H  (b). 


As  may  be  seen  in  the  preceding,  it  is  necessary  to  know  nj©  and 
n®  to  determine  Ix(3^eff  from  R(0). 

As  shown  in  Fig.  1,  the  MFP  for  Rvhv(0)  and  Rhvh(6)  were 
equivalent  (except  for  small  difference  in  amplitude  due  to  the 
Fresnel  factors),  as  expected,  since  all  rays  propagate  with 
refractive  index  no-  The  analysis  using  the  measured  refractive 
indices  gave  lx(3)ivhvAp  ■  Ix(3)Ihvh/<P  with  results  given  in  Table  I.  The 
value  is  smaller  than  lx(3)hso/<P.  as  expected  if  the  components  to 
X<3)  along  the  chain  axis  exceed  those  orthogonal  to  that  axis.  In 
component  notation,  lx(3)lvnv  and  1x(3)1h»h  both  correspond  to 
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FIGURE  3  Maker  Fringe  Patterns  for  a  nematic  solution  of  PBT.  with 
crossed  polars  and  the  incident  polarization  parallel  to  the  director 
arrangement  H-v-V  (a)  and  V-h-H  (b). 

!X(3)!yyyy  with  the  x  and  z  Cartesian  coordinates  along  the  direcior 
and  orthogonal  to  the  sample  plane,  respectively. 

The  depolarized  signals  RHhV(B)  and  Rwh(8)  were  both  nil  (R(0) 
<  0.01),  showing  that  the  signals  Rhvh(B)  and  Rvhv(B)  discussed 
above  are  plane  polarized.  In  the  component  notation  used  above, 
iX'"'l*Hhv  and  lxn)'vvH  both  correspond  to  iX'^'xyyy- 

As  shown  in  Figs.  2  and  3,  the  behavior  is  more  complex  for 
data  on  RwvfO),  Rhhh(B)-  RhvvO)  and  RvhH(9):  the  minima  are 
observed  to  be  nonzero  in  all  cases,  and  for  RhvvIQU  only  one  broad 
peak  is  observed.  Separate  experiments  showed  that  no 

fluorescence  contribution,  caused  by  absorption  at  Xt,^  =  514  nm, 
was  observed  for  a  scan  of  X  from  532  to  800  nm,  permitting 


VM1SOTROPIC  1 ! HKI >  ,  iRl)I  R  M  Si  5  milll-ITY 


TABLE  i  Third-  order  susceonbiiities  tor  PBT 

x-axis  alone  the  ceil  axis 
(the  uirector  tor  nematic  amr- 


1 

Pi 

Ej.ki.  COM  PON  IT 

V  ijj 

i4  ’S13 

ISOTROPIC  SAM  PI  I  .: 

2  ~  \ 

X  (  V  ) 

\(y>  Polarized 

'  M 

x(  y  i 

yi  x  i  Depolarized 

NEMATIC  SAMPLE: 

-i  —  : )  nss 

X 

X  X  X  X  X 

•  o4<  i 

y 

x  y  x  x  x 

x 

y  xvyv 

;  i  | 

<p  is  the  PBT  volume  traction 


elimination  of  this  factor  as  the  cause  of  the  nonzero  minima. 
Although  the  nonzero  minima  ..re  reminiscent  of  die  oehavior 

noted  for  an  isotropic  sample  of  PBT  in  MSA  .it  incident  wavelength 
of  1907  nm.  such  behavior  was  not  observed  with  the  latter  for 

1542  nm  incident  wavelength,  .mu  calculations  using  ::qn  .'>>  ..no 
the  experimental  values  for  a :-»,0  and  show  that  that  is  not  the 
source  of  the  behavior  observed  here.  I'he  nonzero  RuvvPli 
observed  demonstrates  that  x  ’‘  must  ..ave  a  nonzero  c  an 

component  in  the  coordinate  system  used  here,  but  as  discussed 

below,  the  evaluation  of  that  component  may  be  compromised  by 
the  birefringent  nature  of  the  material.  A  manifestation  of  the 
component  off  the  director  axis  is  seen  when  the  intensity  t  ratio 
R  jj  ( 0 . 1- )  * -■  * ict'w rmined  *’  ith  »r!i 1  r •  c ? •: \ n \  '  • h t 

polarized  at  angle  |5  to  the  director  {no  analyzer)  is  studied:  the 
maximum  in  Rp(0,L)  occurs  with  the  ji  *  -  or  -  40  degrees  for  the 
present  monodomain.  This  result  will  be  discussed  subsequently. '• h 
Part  of  the  behavior  is  cieariy  reiaieu  to  me  refractive  index 
for  the  propagating  beams  in  the  birefringent  material.  Thus,  for 
Rwv(8),  the  incident  beam  propagates  with  refractive  index  np  tor 
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Angle .  a 


FIGURE  4  The  parameter  Rc  .0,L)  determined  with  vertical  director, 
vertically  polarized  incident  light  and  analyzer  at  angle  a  to  the 
polarizer  for  sample  thickness  L=385  pm  (top  curve),  and  L=438pm 
(bottom  curve). 


all  0,  but  for  Rhiih(0)«  the  incident  beam  propagates  with  refractive 
index  ne(9),  where 


ne(0) 


jocose  V 

U  °E  , 


'iinifl1'2 

k  n0  J  J 


(8) 


This  accounts  for  part  of  the  difference  in  the  MFP  seen  for  Rvvv(9) 
and  RHhH(9)'.  the  difference  in  the  amplitudes  of  Rvvv(O)  and  Rh>ih(0) 
noted  in  Fig.  3  reflects  variation  in  Kr,  and  is  removed  in  the  final 
analysis.  In  addition,  the  appearance  of  the  MFP  for  Rhvv(6)  and 
RvhH(0)  is  due  in  part  to  the  materials  high  birefringence.  For 
example,  consideration  of  the  values  of  ng^  and  no.3w  shows  that 
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these  are  so  nearly  similar  that  y( 8)  <  it  for  the  available  range  in  9. 
resulting  in  the  broad  peak  observed  for  Rhvv(9)  (Fig.  3a).  The 
increasing  RvhH(9)  noted  with  increasing  0  is  caused  by  the 
decreasing  difference  between  ne>{J>(0)  and  no,3w  with  increasing  9 
(Fig.  3b).  In  addition  to  these  considerations,  a  THG  beam  created 
with  polarization  neither  parallel  nor  orthogonal  to  the  director  will 
propagate  as  two  beams,  one  each  polarized  in  these  directions  in 
the  birefringent  sample,  and  propagating  with  distinct  refractive 
indices  that  may  depend  on  0,  producing  elliptically  polarized 

beam.  For  the  observed  THG  signal,  e.g.,  the  nematic  slab  is  a 

retarder.  This  can  result  in  the  nonzero  minima  noted  in  Figs.  2 
and  3.  A  manifestation  of  this  is  observed  in  Ra(0,L)  determined 
with  vertical  director  and  vertically  polarized  incident  light,  and  an 
analyzer  at  angle  a  to  the  polarizer.  The  maximum  in  Ra(0,L) 

occurs  with  angle  am  that  depends  on  the  sample  thickness  L  as 

shown  in  Fig.  4.  The  difference  in  intensity  between  the  two  cases 
is  due  to  the  dependence  of  the  function  Q(0,L)  on  the  thickness 
(Eqns.  (5)  and  (6)). 

In  order  to  proceed,  apparent  values  of  lx(3)lv»v,  Ix(3),hkh,  Ix(3)|vi>h. 
and  Ix(3),h»v  can  be  deduced  from  Rwv(O),  RHhH(O),  RvhH(O).  and 
Rhvv(O),  respectively,  using  the  relevant  refractive  indices.  These 
are  entered  in  Table  I.  With  neglect  of  the  complications  from  the 
birefringence,  lx(33lvvv  and  lx(3)iHhH»  would  both  correspond  to 
>X(3)|xxxx,  and  lx(3)ivhH  and  Ix(3)Ih»v  would  correspond  to  lx(3)|yxxx- 
Owing  to  the  birefringent  nature  of  the  material,  the  numerical 
values  reported  in  Table  I  for  these  parameters  may  be  suspect, 
but  it  seems  unlikely  that  their  general  ranking  will  be  altered  if 
the  birefringence  is  taken  into  account. 
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/3\ 

Summarizing  the  results,  in  component  notation,  lx  Ixxxx  » 
lx(3)lyyyy  >  lx(3)'yxxx  >  ix(3)Uyyy,  but  lx(3)>xxxx  is  not  the  largest 
projection  of  x(3^-  By  comparison  with  the  data  in  Table  I,  lx  liso/<P 
determined  for  isotropic  solutions  of  PBT  in  MSA  is  given  in  Table  I 
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for  comparison. 10  Of  course,  the  principal  interest  here  i:  not  only 
in  the  lx(3)!ijkl  but  in  the  corresponding  components  lyljjid  of  the 
second  hyperpolarizability  tensor  y.  In  general,  a  relation  between 
X(3)  and  y  may  be  written  as:2-3 

jr(3)(  o;  C0j,C02,C03)  =  f(a>)ffC0|  ;f(a)2)f(CD3)Ny(-0);  (0|. £02,(03)  (  1  1  ) 

where  the  f(tov)  are  local  field  factors  and  N  is  the  number  of 
molecules  per  unit  volume.  An  eventual  goal  is  to  be  able  to 
reliably  estimate  y(-o);  001,0)2,(03)  from  information  on  the  molecular 
structure,  and  to  maximize  the  desired  response  by  manipulation  of 
the  latter.  The  behavior  reported  here  indicates  that  the  maximum 
component  to  y  is  off  the  molecular  axis,  contrary  to  what  one  might 
have  anticipated  A16-17  This  behavior  requires  that  a  relevant 
electronic  transition  in  the  NLO  response  of  the  solvated  PBT  is  at 
an  appreciable  angle  to  the  chain  axis.  This  is  distinctly  different 
from  the  result  of  calculations  on  conjugated  cis  and  trans  polyenes, 
for  which  it  is  predicted  that  the  principal  component  to  y  will  have 
all  electromagnetic  fields  polarized  along  the  molecular  axis.17  With 
PBT,  a  possible  candidate  for  a  vector  associated  with  an  electronic 
transition  off  axis  to  the  rodlike  molecular  axis  is  the  vector  along 
the  sulfur-sulfur  atoms  in  the  repeating  unit  of  the  chain.  This 
possibility  opens  the  question  of  whether  y  for  PBT  is  influenced  by 
the  polymeric  rodlike  character  of  the  molecule,  or  whether  y  would 
be  essentially  the  same  for  a  short  oligomeric  PBT  as  for  a  long- 
chain  molecule,  perhaps  even  a  model  of  the  repeat  unit  (i.e.,  the 
repeat  unit  terminated  by  protons). 
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